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Stochastic  Duels  with  Limited  Ammunition  Supply 

by 

C.  J.  Ancker^  Jr, 


ABSTRACT 

In  a  previous  paper  [l]  the  ground  work  was  laid  for  a  new  theoretical 
model  of  combat  called  the  "stochastic  duel.  "  The  principal  elements  of  the 
model  were  fixed  kill  probabilities  on  each  round  fired^  a  random  time  between 
rounds  fired  and  unlimited  ammunition  supply.  This  paper  extends  the  solution 
to  those  more  realistic  cases  of  limited  ammunition  supply.  The  general 
solution  (as  a  quadrature)  is  obtained  and  it  is  applied  to  several  specific 
examples.  The  special  case  of  flnitej  fixed  ammunition  supply  is  treated  in 


some  detail. 
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Stochastic  Duels  with  Limited  Ammunition  Supply 

by 

C.  J.  Ancker^  Jr. 

In  a  previous  paper  [1]^  the  ground  work  was  laid  for  a  new  theoretical 
model  of  combat  called  the  "stochastic  duel.  "  The  principal  elements  of  the 
model  were  fixed  kill  probabilities  on  each  round  fired^  a  random  time  between 
rounds  fired  and  unlimited  ammunition  supply.  We  now  extend  the  solution  to 
those  more  realistic  cases  in  which  the  ammunition  supply  is  limited.  This 
complicates  matters  somewhat^  but  as  we  shall  see^  it  is  still  possible  to 
obtain  the  general  solution  (as  a  quadrature)  and  to  apply  it  to  specific 
examples. 


ASSUMPTIONS 

The  duel  starts  with  two  contestants  (A  and  B)  both  having  unloaded  weapons 
and  a  limited  ammunition  supply.  A  starts  with  n  rounds.  The  probability  that 
n  is  some  integer  k  is  that  is 

P{n=.k}  -  ,  k.0,1,2,...  (1) 

where 

Jo“k  ■  1  •  «) 

Thus^  the  form  a  discrete  probability  distribution. 

Similarly^  B  starts  with  m  rounds  and 
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P{ni^j)  -  Pj,  j-0,1,2,...  (3) 

and 

jloPj  -  1 

When  the  duel  begins  each  contestant  loads  and  fires  and  continues  this  process 
until  either  one  or  the  other  is  killed  or  they  both  run  out  of  ammunition.  The 
probability  of  a  kill  on  each  round  fired  is  p^  for  A  and  p^  for  B  and  both  are 
fixed.  The  time  between  rounds  is  a  random  variable  described  by  the  probability 
density  function  ^^(0  for  A  and  fg(t)  for  B.  The  solutions  to  this  duel  are  the 
probabilities  that  A  wins  or  B  wins  or  they  tie  (both  run  out  of  ammunition). 

THE  GENERAL  SOLUTION 

The  technique  for  solution  follows  the  pattern  set  in  the  first  paper  [l]. 
Firsts  we  determine  the  distribution  of  the  time  to  a  kill  for  each  contestant. 
These  functions  are  derived  by  considering  the  opponent  to  be  a  passive  target. 
Then  the  probability  that  A  wins  is  simply  the  probability  that  A’s  time  to  a 
kill  is  less  than  B's  time  to  a  kill. 

If  h^(t)  is  the  density  function  of  A's  time  to  a  klll^  then 
P(A  hits  passive  target  B  in  the  interval  (t,t+dt))  « 

\(t)dt  -  a^B(t-«>)dt+a^[p^f^(t)+q^B(t-<»)]dt 

+  . . .  , 
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where  *  denotes  the  convolution  and  6  Is  the  Dirac  Delta  Function.  Thus 


“  Pa 


00  k-1 

k»l  i“o  k«o 


(6) 


oo 

^Al^t)4«(t.»)  , 


k»o 


where  the  superscript  denotes  i  Iterated  convolutions  of  the  function  with 
itself.  Hie  distribution  and  the  complementary  distribution  functions  of  ^^^(0 


are 


and 


u 

"a<')  ■  H*i(0 

^  n 


00  00  00 

°a")  *  Z“k’J  ■  ®Al">  +Z“k4 

^  S _  1 _ 


(7) 


(8) 


k*o 


k«o 


Similar  expressions  apply  to  B  (simply  replace  A  subscripts  by  by  and 

index  k  by  index  j). 

Since  A  cannot  win  If  he  runs  out  of  amnunition  (i.e.  ^  if  his  time  to  a 
kill  is  infinite)  but  he  may  win  if  B  runs  out  (i.  e. ,  A  continues  to  fire  at  B 
until  out  of  anmunition)  ^  the  probability  that  A  wins  is 


P(A)  - J Gg(t)h^^(t)dt 

00  00  00 

-y'Ggi(t)hAi(t)dt  +^Pjqjy'h^^(t)dt  , 
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or  from  (6) 


00  00  00 

P(A)  -  r G3^(t)h^^(t)dt  +  [l  •  (9) 

°  j-o  k-o 

From  Parseval’s  formula  [2]^  equation  (9)  may  be  written  as 
^  00  00  00  00 

.00  °  j-0  k-0 

where  characteristic  function*  of  and  is  defined  as 


00  k-1 


00  k-1 


k-1  £-0 

Pa  Z^kl^fcpf^u) 


(11) 


k*l  i*o 


00 

k*o 


and  where  Cp^(u)  characteristic  function  of  f^(t)  or 

00 


(12) 


%^l(u)  la  not  a  proper  characteristic  function  since  ^  proper 

density  function.  However^  the  definitions  and  manipulations  given  are  valid. 
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In  the  derivation  of  (11)  the  convolution  property  of  characteristic  functions 
has  been  used.  ^ 

CO 

The  inner  integral  of  (10)  integrates  by  parts  to  ■^|^<I>gj^(u) »!+  ^  j 
\diere  Is  defined  as  in  (11)  and  is  ^  ^ 


*BiW  • 

j-0 


(13) 


and  cp«(u)  is  the  characteristic  function  of  f«(t).  Thus 


00  00 

p<A) .  [i 


j»o  k«o 

+  00 


s/  *A1 

-  00 


(-U) 


$  ,(u)  -  1  H-Vp 


jwQ 


1  jl  du 

’j’bJ  -■ 


(14) 


It  may  be  shown  that  functions  such  as  <I»(u)  and  cp(u)  have  the  following  properties: 

(1)  no  singularities  in  the  upper  half  of  the  complex  plane  and  (2)  diminish  like 

~  along  a  semicircular  contour  of  radius  R  in  the  upper  half  of  the  complex  plane. 

R 

<t(-u)  and  <p(-u)  have  similar  properties  in  the  lower  half  of  the  complex  plane. 
Sufficient  conditions  for  the  second  property  are  that  f(t)  be  a  differentiable 
function  of  bounded  variation  and  lim  f(t)  «  0.  From  (11)  and  (13)  we  also  note 


that 

00 

k*o 

GO 

j«o 


\ 


/ 


(15) 
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In  (14)  the  term  - 


r  ^ 

2jci  J  X 


•u)du 


r 

-  /  - -  may  be  integrated  by  noting  from 

>  00 

the  foregoing  that  in  the  complex  domain  the  integrand  has  no  poles  in  the  lower 
half  plane,  the  integral  vanishes  around  the  contour  C  as  R-w  (see  Figure  1), 

and  there  is  only  a  simple  pole  of  residue  - 


1»o 

2ici 


G'Iva) 

k«o 


at  the 


origin  (see  equation  (15)).  Thus  letting  p-+0  (and  noting  that  we  take  - 


Figure  1. 


times  the  residue  because  the  pole  is  outside  the  small  semicircle),  the  desired 


integral  is  ^  ^1 

k«o  j«o 


00  00  +00 

P(A)  ■  i[l  -^Va]  ♦Ipj’b]  *■  «Al( 

k*0  j»0  -  00 


-0)4181  (“) 


du 


(16) 
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The  Integral  in  (16)  is  the  Cauchy  principal  value.  This  may  be  further  simplified 

by  noting  that  the  integrand  in  (16)  has  a  simple  pole  at  u"0  which  contributes 
00  00 

to  the  integral  (path  of  integration  as  in  Figure 

k»o  j«o 

1)  so  (16)  may  be  written 

00  00 

P(A)  [i  f .  an 


k*o 


where  J  means  integration  from  <-«  to  -H»  along  the  real  axis  excluding  the  simple 
pole  at  the  origin.  Further^  it  implies  that  the  path  of  integration  is  as  in 
Figure  1  (where  R->«>),  and  that  the  integrand  is  analytic  everywhere  except  at 

V 

singularities  and  that  the  integrand  vanishes  at  least  like  7^  on  C  (where  K 

R 

is  some  constant)^  and  that  p  is  finite  but  less  than  the  distance  to  the  nearest 
singularity  in  the  lower  half-plane.  The  direction  of  integration  is  in  the 
negative  sense. 

A  still  different  form  may  be  obtained  by  integrating  around  the  path  of 
Figure  2.  In  this  case  the  pole  at  the  origin  contributes 


Figure  2. 
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C 


00  00 

k«o  j»o 


00 

P(A)  -  1  ^ 

k»o 


(18) 


where  /  means  integration  from  -oo  to  4c»  along  the  real  axis  excluding  the  simple 

Ju 

pole  at  the  origin.  Corresponding  restrictions  apply  to  this  path  of  integration 
in  the  upper  half -plane. 

In  the  event  the  integrand  is  not  analytic  or  does  not  behave  properly  on 
P(A)  must  be  evaluated  from  (16)^  rather  than  (17)  or  (18). 

Similarly^  the  probability  that  B  wins  is^ 

00  00  +  <»* 

m  -  i[i  [i  f 

j«o  k-0  -  «> 

00  00 

■  I Va  /„  »ai<-“)‘bi<“>  V  («) 

k»o  J»o 

00 

"  ^  -1^)4  •  2^  /.  ♦ai(-“>*B1<'‘>  It  • 

j-o 


In  addition^  the  probability  of  a  draw  (the  probability  both  run  out  of 
ammunition  before  a  kill)  is 


00  00 


KAB)  .  , 


(20) 


k*o  j«o 
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and,  of  course,  P(A)  +  P(B)  +  P(AB)  ■  1. 

Thus,  (16),  (17),  (18),  (19),  and  (20)  along  with  (11),  (12)  and  (13)  con¬ 
stitute  the  general  solution  to  the  stochastic  duel  with  limited  annunition  supply. 

These  results  specialize  to  the  unlimited  ammunition  supply  case  of  the 
previous  paper  [l]  by  letting  -  Pj  -  0  for  all  finite  k,J  and  -  1 

for  k, j  infinite. 

SOME  EXAMPLES 

In  this  section  the  general  solution  is  applied  to  certain  examples.  In 
each  example  there  is  a  specific  distribution  for  the  number  of  rounds  A  starts 
with,  the  number  of  rounds  B  starts  with,  and  specific  density  functions  for 
their  times  between  rounds.  These  are  Indicated  in  the  heading  for  each  example. 

In  the  following,  only  P(A)  and  P(AB)  are  given  since  P(B)  -  1  -  P(A)  -  P(AB). 
Example  1.  Geometric -Geometric-Negative  Exponential 

Let 

-  (l-a)c^,  pj  -  (l-p)pJ 

"  V  >  fB(^>  -  V 

where  the  parameters  r^  and  r^  are  average  rates  of  fire.  Then, 

■  7^  ■  tB*")  ■  7^ 


and  from  (11)  and  (13) 
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PPb"b 


®A1^''‘^  "  (l-aq^)r^+iu  '  •bI^'*^  "  (l-pq3)rg-iu 


Thus,  from  (17) 


p.M  ^  rZ—^Af^'N  /  PPb’^B  N  ^ 

“  (l-pq^Xl-Oq J  2«i  jA(l-aq.)r  +iuy''^(l-pq  )r  -iu;  u  ’ 


PPp>^r 


'A' ‘A" 


'B'  B 


The  Integrand  has  one  singularity  in  the  lower  half-plane  which  is  a  simple  pole 

1  '^PaPb’^a 

.t  u  .  -ici-pq,)r,  with  residue  -  —  (I.^-)T(i-S?;)r^l.(l-p,^)r,-]  ' 
applying  the  residue  theorem 

OPa  r  (l-QqA)rA+(l-P)rB 


P(A) 


°^A  U-QqAjrA+ii-p;rB  . 

“  l-oq  [(l-Qiq  )r^+(l-pq  )r  ' 


(22) 


From  (20), 


p(ab)  - 

(l.aqj(l-pqj  • 


(23) 


In  the  following  examples  only  the  assumptions  and  the  solutions  are  given  since 
the  procedure  for  derivation  Is  the  same  as  In  Example  1. 

Example  2.  Poisson-Geometric -Negative  Exponential 


°k  -  ^ 


and 
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P(A) 


r,  .FPaV^^-pV’^bI  1 

J i 

(l-^qg)  LPaV^^-PV’^B^ 


(l-.-'*A) 


(24) 


and 


P(AB) 


(lzP)e'°^A 
l-Pqfi  • 


(25) 


The  Geometric-Poisson-Negative  Exponential  case  may  of  course  be  obtained 
by  simply  reversing  the  roles  of  A  and  that  is  exchanging  A  and  B  subscripts, 
and  replacing  a  with  ^  and  vice  versa. 

Example  3.  Binomial-Geometric -Negative  Exponential 


/  1  Y  /  k  4 

^  "  \i45y  \  k"0,l,...,n;  pj  ■  (1“P)P 

fA(t)  -  r^e  ,  fgCt)  -  r^e 


P(A)  - 


PPaPb’^a 


-  fl  -  ^  [l  ,  T 1 

>  t  L  V<^-P‘Ib>^bJ  j 


(26) 


and 


.  {^')  -UiEL 

\l+a  J  (l-pqj) 


P(AB) 


(27) 
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t 


Again  the  Geometric -Binomial-Negative  Exponential  case  may  be  obtained  by 
reversing  the  roles  of  A  and  as  stated  in  Example  2,  and  replacing  n  with  m. 
Example  4.  Geometrlc-Geometric-(y ,2) 


0^  -  (l-a)Q?',  pj  -  (l-p)pj  , 


,  -2r  t  ,  -2r  t 

f^(t)  -  4r^e  ,  fg(t)  -  4rge 


p(A)  -  r _ (l-oqA)r;-(l-pqB)r^44rB(r^+rB) _ > 

H(l-aq.)r?-(l-pq_)r?]^+4r.r_(r.+r_)[(l-aq.)r.+(l-pq_)r_]J 


■‘A'‘A 


'B'  B'  '  ‘A  B'  A  B' 


'‘A' ‘A" 


»B'  B' 


Q!(1-P)Pa 

+  (l-oq^xi-pq^y  ^ 


(28) 


and 


i 


P(AB) 


d-caa-B] 

(I-Ctqjil-Pqg) 


(29) 


The  Important  special  case  of  a  fixed  limitation  on  ammunition  supply  will 
now  be  ^.onsidered  in  more  detail. 


THE  CASE  OF  FINITE,  FIXED  AMMUNITION  SUPPLY 

Suppose  that  A  always  starts  with  exactly  n  rounds  and  B  with  exactly  m 
rounds.  This  might  very  nearly  be  the  case  where  the  policy  is  to  disengage 
and  resupply  immediately  after  action  and  where  there  is  a  physical  limitation 
on  the  number  of  rounds  which  can  be  carried  such  as  in  the  Interior  of  a  tank 
or  an  interceptor. 

We  may  go  at  once  to  the  general  solution  for  this  case  by  letting 
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0^-1,  k  -  n 

*  0  all  other  k 


Pj  -  1,  j  -  m 


■  0  all  other  j. 

Thus,  from  (11),  (12),  (13),  (16),  (17),  (18),  (19),  and  (20) 

+  “ 

P(A)  f 


and 


■  ‘-’I  *  2h/„‘ai<-“)»bi<">  ^  ' 

+  “ 

P(B)  -  ®^i(-u)®3i(u)  ^ 


'  ife  /u®Al(-'‘>®Bl('‘>  f 
■  2*1  X'^Al^'^^^Bl^'*^  T 


(30) 


(31) 


(32) 


and 


P(AB)  -  , 


(33) 


where 
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■  TvIfT  L‘  ■  (''a»a<“V  J  ' 


and 


pr<pb(“)  r  /  Yi 

"  l-qgCPg(u)  ■  (Vb^'‘7  J  • 


(34) 


(35) 


I’f 

Jl 


and  the  cp's  and  all  have  the  same  meanings  and  restrictions  as  before. 

The  use  of  the  general  solution  (equations  (31)  through  (35))  will  now  be 
illustrated  by  an  example.  Assume  that  firing  times  are  negative  exponential^ 
thus 

^A<‘>  "  V  ^  “  V  ® 


From  (31) 


m,,  n^  .  Pa^b'^A^’B 
P(A)  .  q^d-q^)  + 


II— 


/  •  J  \  ^  n  n 

(r^+iu)  -q^r^ 


A  m  mm 


(Pa'^A+^«)  (’^A+i^) 


«  ,  V  II 

r  "Vb 


T  • 


The  integrand  of  (36)  has  only  one  pole  in  the  lower  half-plane.  This  is 

located  at  m’th  order  pole.  The  residue  at  this  pole  is 

1  d®*^  m 

given  by  Res  *  (^-1)7 — jJJTi 

'  '  du  ^o 

integrand  of  (36).  From  this 
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8(u) 


2jti 


,n  _n  n 


PaPbVb  ^j-  1  r 


_i  'I 

L(p,r„-iu)(.l)“uJ  J  ■ 


■(Pa>^A'''^“^^’^A'''^“^  -■  ^(Pb^b-IuX-I) 
The  (m-l)8t  derivative  of  this  reduces  to 


,m-l 


du 


u»u^ 


vin  d'  rw  ^  i 

— „  —  (iqsrB)  -^[h(u)] 


du 


u«u 


where 


'  (r*+iu)"]  ■ 


(r^+iu)* 

Expanding  h(u)  by  partial  fractions. 


h(u) 


n  n 

ri._v^ir._ 

L  fr.+iu^^^J  L  Pa^, 


(r^+iu)*‘J  L  Pa^A<Pa’^A+Pb^B><PaV^'^> 

+ _ i _  ,  1 

Pb’^b^PaVPb’^b^  Cpb’^b-^'^>  PaPbVb'*. 


The  (m-l)8t  derivative  of  (40)  may  now  be  taken  by  the  product  formula 
is  Inserted,  and  the  whole  expression  substituted  into  (38).  Finally, 
to  the  residue  formula,  we  have  after  much  calculation 


(37) 


(38) 


(39) 


m 


Then  u 

o 

going  back 


23  April  1963 


18 


SP- 10 17/00 1/00 


This  expression  may  be  put  In  closed  form  by  using  the  known  [3]  relation 


k-o 


(l-x)‘”[l-I  (m,n)]  , 


(42) 


where  l^(m,n)  Is  the  well-tabled  [4]  Incomplete  Beta-Function  Ratio  defined  by 


B  (m,n) 


-  r(m)r(n) 


A 

r(nHii)  r  m-1,.  .n-1 

'(m)r(n)  JJ  • 


The  Identity 


(43) 


is  also  used.  When  (42)  and  (43)  are  applied  to  (41)  the  final  result  is  obtained 
after  some  manipulation. 


(44) 


where 


^  - 


(45) 


To  complete  the  example^  from  (33), 


P(AB)  -  q"qj  . 


(46) 


Two  other  cases  are  obtained  rather  easily  from  these  results.  If  A  has  a 
fixed  number  of  rounds^  and  B  has  an  infinite  number  of  rounds^  then  for  negative 
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exponential  firing  times,  the  solution  is  obtained  by  letting  m-**  in  (44),  (45) 
and  (46). 


P(A) 


pa’^a  r 

1  ( 

'Va  Y 

'1 

Pa’^A'^Pb’^B  - 

^Pb’^b+V 

(47) 


and 


P(AB)  -  0 


(48) 


Similarly,  if  A  has  an  infinite  number  of  rounds  and  B  is  restricted,  let 


n-^. 


Pa^a  ,  /  Pb^b  /  Vb  Y 
“  Pa’^a+Pb’^b  Wa+PbV  WVV 


(49) 


and 


P(AB)  -  0  . 


(50) 


In  the  preceding  we  have  used  lim  I  (m,n)  ■  0  and  1^  I  (m,n)  ■  1 
^  m-MB  *  n"**®  * 


DISCUSSION  AND  CONCLUSION 

The  general  solution  to  the  stochastic  duel  with  limited  ammunition  supply 
has  been  obtained.  Several  examples  with  specific  distributions  have  been  worked 
out.  The  general  solution  for  the  special  case  of  fixed,  finite  ammunition  supply 
has  also  been  obtained  and  several  examples  worked  out. 

pggults  given  here  can  be  readily  adapted  to  fit  various  initial  con* 
dltlons  as  described  in  the  previous  paper  [l].  These  Include  the  "Classical 
Duel"  where  each  contestant  starts  with  loaded  weapons  and  the  "tactical  equity" 
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duel  in  which  one-half  the  time  A  sights  B  first  and  fires  one  round  which  alerts 
B  and  then  the  duel  proceeds  as  a  fundamental  duel.  The  other  half  of  the  time 
B  fires  first. 

A  limited  ammunition  supply  may  drastically  affect  the  outcome  of  the  duel. 
The  precise  effect  depends  on  the  exact  conditions  of  the  duel.  Some  insight  into 
the  nature  of  this  effect  is  clearly  shown  from  an  examination  of  equations  (47) 
and  (49).  These  are  duels  in  which  firing  times  are  exponential  and  only  one 
side  has  a  fixed^  limited  ammunition  supply.  We  now  compare  these  results  with 
the  same  duel  in  which  both  sides  have  unlimited  supplies  whose  outcome  [l]  is 


P(A)u 


Pa’^a 


PaVPb’^b 


(51) 


In  (47),  A  has  a  limit  of  n  rounds.  The  effect  of  this  limitation  is  given 


(52) 


This  expression  is  plotted  in  Figure  3.  A's  chances  of  winning  are  always  degraded 
for  n  <  00,  If  p^r^  is  much  larger  than  r^,  then  x  is  nearly  zero  and  P(A)/P(A)y 
is  nearly  1  and  the  limitation  on  ammunition  has  little  effect.  This  is  explained 
by  the  fact  that  P(A)  is  very  small  in  this  case  anyway.  If  p  r.  is  of  the  same 
order  as  r^  or  is  much  less  than  r^  then  x  essentially  depends  on  (A’s  miss 
probability).  If  is  small  then  again  the  effect  is  small  but  if  q^  is  fairly 
large  then  the  ammunition  limitation  is  important  and  A’s  chances  of  winning  are 
degraded  seriously  by  small  n. 
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P(A) 

P(A)y 


0 


0.5 


1.0 


X 


Figure  3.  The  Effect  on  A  of  Ammunition  Limitation  on  A. 


For  the  duel  where  B’s  supply  is  limited  (equation  (49))  we  have 

1-P,  lm 


P(A)u 


1  + 


Vb 

Pa’^a 


1+P 


B 


aIa 


(53) 


This  expression  is  plotted  in  Figure  4.  A's  chances  of  winning  are,  of  course, 
always  Improved  for  m  <  ».  However,  the  major  effect  occurs  for  small  p^  and 


this  is  particularly  reinforced  for  p^r^  <  Pg^g. 
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APPENDIX 

The  general  solution  for  a  somewhat  different  duel  is  readily  obtained  from 
the  preceding  work. 

In  this  case^  the  duel  concludes  as  soon  as  one  contestant  is  killed  or  as 
soon  as  one  contestant  runs  out  of  ammunition.  The  last  feature  is  new.  This 


might  be  the  case  where  the  contestant  who  runs  out  takes  immediate  cover.  The 

analysis  proceeds  as  before  but  account  is  taken  of  the  fact  that  no  kills  are 

possible  when  one  contestant  has  no  ammunition.  The  modified  general  solution  is 
00 

P(A)  - J G3^(t)h^^(t)dt 

00  00  ^00 

k«o  j»o  -  « 

(54) 

“  2^1  ^ 

00  00 

■(‘  -Ivi)  0  /ai<-")»biW  r  > 

k-O  j-0 

and  «  oo  +00 

P(B)  .  i(l  (l  ^  f  f 

k«0  j«0  -  00 

■  -  Sr  /„‘ai(-“>«B1<“)  T  <55) 

00  00 


00  00 


p<AB)  • 

k«o  j«o  k»o  j«o 


The  ®’s  are  the  same  as  before.  All  of  the  examples  worked  out  in  the  preceding 
sections  are  easily  modified  to  fit  this  model. 
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